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ABSTRACT 

The lensing properties of the Plummer model with a central point mass and external shear are 
derived, including the image multiplicities, critical curves and caustics. This provides a simple 
model for a flattened galaxy with a central supermassive black hole. For the Plummer model 
with black hole, the maximum number of images is 4, provided the black hole mass is less 
than an upper bound which is calculated analytically. This introduces a method to constrain 
black hole masses by counting images, thus applicable at cosmological distance. With shear, 
the maximum number of images is 6 and we illustrate the occurrence of an astroid caustic and 
two metamorphoses. 
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1 INTRODUCTION 

Albert Einstein, in his 1911 paper On the Influence of Gravitation 
on the Propagation of Light, applied the principle of equivalence 
to the deflection of light rays by the sun. This marks the beginning 
of our modern understanding of gravitational lensing. After the full 
prediction from General Relativity had been confirmed by Edding- 
ton's observations of the 1919 solar eclipse, it was only in 1936 
that Einstein considered, rather pessimistically, the deflection due 
to stars. Zwicky's (1937) insight one year later that galaxies were 
much better candidates was finally corroborated by the discovery 
of the first doubly-lensed quasar in 1979 by Walsh, Cars well and 
Weymann. 

Nowadays, over a 100 multiply imaged gravitational lens sys- 
tems are known. In many of these, a background quasar or ra- 
dio galaxy is lensed into doublets or quadruplets by an early-type 
galaxy. Provided the density distribution at the centre of the lens- 
ing galaxy is weakly singular (that is, no worse than p ~ 1/r), 
then gravitational lensing yields an odd number of images, gen- 
erally three or five. If the singularity is stronger than this, then 
there is an even number of images, as the central image is absent 
lEvans & Wilkinsonll 9981. Most early-type galaxies are believed 
to harbour supermassive black holes in their centres. The existence 
and the magnification of the central image in gravitational lens sys- 
tems are very sensitive to the mass distribution in the centre and the 
presence of a black hole. lKeetonl <2003l) has shown that the likely 
range of the flux of any central image is below the range of cur- 
rent instruments, but within the grasp of forthcoming, high sensi- 
tivity VLBI arrays. Thus, it is becoming possible to probe the mass 
distribution in the central tens of parsecs of the lensing galaxy. It 
is therefore worthwhile to de velop some s imple models of lensing 
by galaxies with black holes. iMao et al"l 1200 ll) already studied a 
galaxy model akin to the cored isothermal sphere and showed that 



a central black hole - provided it is not too massive - can create an 
additional faint central image. 

The present paper is a mainly analytical study of the Plummer 
lens with a central black hole and in the presence of external shear. 
In §2, we briefly review the properties of the spherical Plummer 
lens and show that its images obey a magnification invariant. This 
model is extended in §3 to include a central point mass to which 
external shear is added in 54. 



2 THE PLUMMER MODEL 

Plummer's (1911) model is the simplest self-consistent solution 
of the collisionless Boltzmann equation with finite total mass. Al- 
though the density in early-type galaxies decreases less steeply than 
p ~ r~ 5 at large radii, the Plummer model still gives a good de- 
scription of the inner parts, which control the lensing properties. 

For an isotropic stellar system, the distribution function / de- 
pends on the (negative) total particle energy E only. The Plummer 
model is then a polytrope defined by the distribution function 



f{E) oc (-E) 



7/2 



(1) 



For a static system, energy conservation E — <f> + v 2 /2 holds and, 
as a consequence, Poisson's equation for a self-consistent, self- 
gravitating, spherical and isotropic system 



LA. 

r 2 dr 



16tt G / f(E)v 2 dv 



(2) 



implies the mass density (e.g. Binney & Tremaine 1994, p. 225) 

,2\ ~ 5 / 2 

ro 



p = po\l + 



(3) 
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and total mass M = 47rpo?'o/3. 
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Figure 1. Plummer lens with central point mass, kq = 2, ao = 0.02. Left: critical circles in L, right: caustic circles in S. Number labels indicate number of 
images. 



2.1 Lens Equation 

The angular diameter distances from the observer to the lens plane 
L, the source plane S and from L to S be denoted by Dlo,Dso 
and Dsl, respectively. We use the conventional dimensionless co- 
ordinates x in L and y in S (see e.g.. lSchneideretal. Ill999l p. 
245). The relative surface density ftp is now obtained by projecting 
onto L, 



ftp 



fto 



(-y 



(4) 



where x = |x| such that the origin coincides with the lens centre, 
x = r /I?LO and 



M 



fto = 



Scrit 



c'Dso 



4ttGZ)sl.Dlo 



The deflection potential ^ is 



fto^o , / , fx 
*p = *o + In ( 1 + ( — 

x 



The lens equation is then 
y = x - V*p = x 



K x 



l + ( — ) 2 ' 



(5) 



(6) 



(7) 



The image and caustic properties of the Plummer lens are briefly 
summarized below, for comparison with our results in §3 and §4 



2.2 Image and Caustic Properties 

The determinant of the Jacobian of the Plummer lens map is 

2ft (^) 2 



det J : 



fto 



+ 1- 



l+(-2L)2 I (1 + (^)2V 



•(8) 



The origins of L and S are critical and caustic points, respectively. 
For a point source at the origin y = 0, the lens equation implies a 
circular image, the Einstein ring, with radius x — xoV^o — 1 and 
infinite magnification. Otherwise, without loss of generality, let the 
source position be y = (yi,0), whence x-z = and 

x\ - y-ixi + (1 - K )xlxi - x 2 y± = 0. (9) 

Thus, there are one or three images, depending on the source po- 
sition yi , The images lie on a straight line through the lens centre. 



For — > co, J — > I (the identity matrix), so there is always one 
image at large yi , as expected. 

The critical curves can be obtained from {8j with det J — 
resulting in 



aJcrit.l = iovko 



£crit,2 




There are two critical circles in L if ko > 1. Applying the lens 
equation then shows that the caustic corresponding to x C rit,i is the 
origin point in S, whereas the other caustic is circular, as expected 
by symmetry. Hence three images can occur by crossing this caustic 
inwards. However, if fto ^ 1 then only a critical and caustic point 
exist and the number of images is one for all y / 0. 

2.3 Magnification Invariant 

Witt & Mao's (1995) paper introduced the notion of a lensing in- 
variant. They considered a binary lens of two point masses. Their 
main result is that the sum of the signed magnifications of the max- 
imum number of five images is unity, namely 



5 

E 

i=i 



= i. 



(10) 



Here, \ii is the absolute magnification of the image, while pi is 
the parity. This result holds good irrespective of the position of the 
source, provided it remains within the central ca ustic giving rise to 
the maximum numb er of images. Subsequently, RhieJ 119971) and 
IWitt & Mao I J2000) found further examples of lens models with 
invariants. A gene ral theory of such invariants using sing le variable 
complex analysis iHunter & Evans Il200ll : lEvans & Hunter ll2002t 
lAn & Evansll2006l) was developed. This approach is used here to 
investigate the general circularly symmetric lens and, in particular, 
the Plummer model. 



2.3.1 Circular Symmetric Lenses 

The standard complexification is introduced (e.g. Petters, Levine 
& Wambsganss 2001, chap. 15.1). Let L c , S c be the complex lens 
and source planes, respectively, and define the complex coordinates 



xi + ix2 € L c 



C = yi + W2 s s c 



(ii) 
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Let the complex conjugate be denoted by a bar. The complex lens 
equation is rendered 



C = z- zf (zz) , 
where the deflection factor / is 
,d*(a) 



/(-) = 2- 



c/.s 



(12) 



(13) 



(15) 



Note that, by definition, / is a real-valued function so that / = /. 
The conjugate of the lens equation is therefore 

( = z-zf(zz) (14) 

and so 

-f 

The geometrical meaning of 1151 is plain: due to circular symmetry, 
we can without loss of generality orientate the real axes in S c , L c 
so that the source position f = f is real and hence, by 1151 . the 
image positions as well: lens centre, point source and all images are 
collinear. Also, 1151 allows us to eliminate z from the lens equation 
1 1 21 so that standard, one-dimensional complex analysis may be 
applied. It can be recast thus to give the imaging equation 

J(*,£C)=0, (16) 
defined so that image positions z< correspond to the roots of 



Z(z,C,C) =z-zf 



c 



(17) 



As the lens model / may have isolated poles in the complex plane, 
the imaging equation is meromorphic. 

The signed flux magnification of the image at z% is given by 
(e.g.. |Petters etal.l200lL p. 85) 



(18) 



det J(zi) ' 

where the Jacobian J and its inverse J -1 in complex form can be 



i)z 



dz 



dz I _ 

5C k 



dz I 



along a contour C enclosing all A max roots m . Assuming that the 
roots of the imaging equations are simple, the integrand F has a 
simple pole at each root of / and hence residues 



s(F;Zi) = 



(23) 



by Taylor expansion of I about Zi, as required. Now from the 
residue theorem, we find 



N 

E 

1=1 



N 

E 

<=i 



res(F- Zi 



1 

2^ri 



{C + r)i 



dz. 



(24) 



Furthermore, two assumptions are made which both restrict the 
functions /. First, no branch cuts to ensure singlevaluedness of I 
are to occur, such that the complex plane is analytic everywhere 
outside C. Secondly, it is assumed that the deflection angle becomes 
negligible at large distance from the lens centre. This amounts to 
V 1 !' — > as x — > oo, hence zf — > and therefore I — > z as 
\z\ — > oo. Then by Cauchy's theorem, C can be distorted to a circle 
at infinity Coo 



JV m „ 

E ^ 

i = l 



1 



« + !)! 



dz 



1 



dz 

z 



(25) 



This result is a magnification invariant for the maximum number of 
roots iVmax of the imaging equation. 

Although each image corresponds to a complex root of the 
imaging equation, there may exist roots which do not correspond to 
physical images and which are therefore called spurious roots. To 
understand their impact on 1251 . consider two of their properties. 
First, by choosing without loss of generality £ = ( in the imaging 
equation 1161 . it is seen that both z and z are solutions. So, all 
spurious roots occur as complex conjugate pairs. Secondly, from 
(TJj) and J2j}, we see that the magnifications of a pair of conjugate 
spurious roots are also complex conjugate. This implies that the 
total contribution from spurious roots to the signed magnification 
sum is always real. 



<■).: 



dz\ z 



dz\_ 



dzl 



Requiring that J J 1 = I, it follows that 
detJ=-g^. 

OZ 

5ck 

By standard calculus, we find that the sum of the signed magnifica- 
tions of the N images is 



(19) 



E^ = E 



det J(zi) 



E 



9/1 



ax I ac 

~1 dz | <* C 



(20) 



By using i 1 21 in the denominator and i 1 61 in the numerator of I20i . 
the following simplification is obtained 



E w= E 



"> dz 



(21) 



We now seek a complex-valued function with residues of the form 
J21L Consider the integral of 



F(z,C,0 



(22) 



2.3.2 The Particular Case of the Plummer Model 
The Plummer lens equation can be rewritten as 

K Z 

c = z -TTTz- 

From the definition of the deflection factor /, it follows that 

Ko 



f(z) 



1 + 



C 



(26) 



(27) 



Hence zf — *■ as \z\ — > <x and / involves no roots of z and 
hence no branch cut. The imaging equation is of third degree in 
z and, since the maximum number of images is three, there are 
no spurious roots. Thus, the Plummer model has a magnification 
invariant 



(28) 



In other words, the sum of the signed magnifications [lipt of the 
images is an invariant of the model. This result appears not to have 
been realized before. 
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3 THE PLUMMER MODEL WITH BLACK HOLE 
3.1 Lens Equation 

The standard Plummer lens is now modified by introducing a point 
mass m as a model for a black hole at the centre. Using dimension- 
less coordinates in L as before, its relative surface density is 



kbh(w) 



m e(2) 



<5 (2, (w), 



(29) 



whence the total deflection potential is 



*(x) = * P (x) + i / K B H(w)ln|w-x|d 2 u> (30) 

= * P (x)+alna;, (31) 

where a = m/7r7?L0 2 Scrit- The lens equation follows immedi- 
ately, 



y = x 



kqx ax 



l + (— ) 2 x 2 



(32) 



using 



3.2 Image and Caustic Properties 

This lens equation shows that the point mass introduces an addi- 
tional image at x — > —ay/y 2 as y — > oo: the minimum number 
of images is two. Moreover, eq 1321 is of fourth degree in x, so a 
maximum of four images is expected. Some further insight can be 
gained analytically. With the substitution 



Xo 



= £ cos ( 



a-2 

Xo 



a m 

x 



7rrgE cri t 



(33) 



where £ ^ 0, the Jacobian can be constructed from the partial 
derivatives, 



dyi 
dxi 

dyx 

dX2 

dm 

3X2 



f2 , 



2/t £ cos 2 6> ao(2cos 2 <9- 1) 

(i+c 2 ) 2 + e 



2kq£ cos 9 sin 9 2ao cos 9 sin 9 



= 1 - 



(l + C 2 ) 2 



2ko£ 2 sin 2 6 ao(2sin 2 



9y2 
dx i 



i + £ 2 1 (l + C 2 ) 2 1 e 2 

The radii of the two critical circles are given by det J = 0, which 
yields: 

K do 



1 1 + f 2 f 2 

- 1 - ' Scrit,l Scrit,l 

so that 



= 



and 



\ 



«o 



«0 + do — 1 



(/to + o ) 



+ a 



2 K o£crit,2 . OO 
i T~t , To To T 



1 + £crit,2 (1 + Crit,2) 2 fcrit,: 



0. 



(34) 



(35) 



(36) 



The critical circle I35i has precisely one solution for all positive 
/to, Oo. Applying the lens equation shows that it maps to a caustic 
point at the origin of S. If ko ^ L there is hence one critical circle 
in L at £ cr it,i and one caustic point in S: the number of images is 
always two for y 7^ 0. If no > 1, £crit,2 introduces two further 
critical circles in L and two caustics in S, in contrast to one for 
the Plummer model. Although 1361 is a third degree polynomial in 
£crit,2 and could in principle be solved analytically, it is, however, 



more instructive to display the critical circles and caustics numeri- 
cally. FigureQillustrates this case with three critical circles and two 
caustic circles, giving rise to an annular domain with four images 
inside and two outside of it. 



3.3 Maximum Black Hole Mass 

If the point mass ao increases while the Plummer mass /to is kept 
fixed, the former will eventually dominate the latter, for some ao > 
oo,max- In consequence, two images will occur for all y 7^ 0, as 
is plain from the lens equation for negligible /to. The annulus in S 
resulting in four images will hence disappear, as shown in Figure 

HI 

Let us take without loss of generality y = (yi , 0) =>■ X2 = 0, 
then 1321 requires for the presence of four images Ko > 1 and four 
intersections of 

f vi go6 , o 

Si = i 1 C 2 + T- (3') 

x 1 + ?i £1 

where £1 = xi/xo- At the threshold £1 = £0, 00 = 00, max, 
the graphs of left-hand side and right-hand side match in first and 
second derivatives, 



Ko(l-£o) o , ; 



= 



2/to£o(£o — 3) 2ao,n 



(i + £o 2 ) 3 



So 



(38) 
(39) 



Eliminating ao. max , they yield a third degree polynomial in £0 
which becomes, after the standard substitution q = £q + 1, 



q + 3/tog — 4/to = 0. 



(40) 



Since its cubic determinant is positive, there is one real solution 
which can be obtained from Cardano's formula, 

q = (2/t + ko\/^T+4) 1/3 - (kqVko+4 - 2k ) 1/3 (41) 
to give, using <39> and J40I . 

\2/ 



ao 



(o-l) 2 (4- 



(42) 



4-3g 

In physical units, the maximum mass of the point mass for four 
images relative to the Plummer model is 

^max Oq max 



M 



«'■() 



(43) 



Now consider this as a model of a galaxy with a supermassive 
black hole. Since many parameters are free, let us take, as a typ- 
ical example, the total galaxy mass to be M = 10 12 Mq and 
scale length ro = 10 kpc. Let the lens and source be at red- 
shifts 2 = 0.5 and z — 3, respectively, in an Einstein-de Sitter 
universe with H = 70 kms -1 Mpc -1 . This is a sufficient ap- 
proximation to the Concordance Cosmology for the purpose of es- 
timating parameters here, especially because the source distance at 
higher redshift almost cancels out. A well-known expression for an- 
gular diameter distances (e.g., Peacock 2003, p. 103), then yields 
Dso, -Dsl, Dlq. Hence from the definitions we have: 



AGM D SL D LO 
ko = —5 =; = 1-28 > 1 

c^ro Dsoro 



(44) 



so that four images can occur, provided the mass of the central 
black hole is less than 



max = 4.2 x 10~ 3 



M 



= 4.2 x 10 9 M 



(45) 
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Figure 2. Plummer lens with central point mass, reo = 2. Left: radii of critical circles in L as function of point mass ao, right: radii of caustic circles in S as 
function of ao- Number labels indicate number of images. Note the disappearance of the caustic annulus with four images. 



from I43i . In fact, this result is of the same order of magnitude as 
the masses inferred from observations: the presently surmised me- 
dian mass of supermassive black holes is (1.4 ± 0.4) x 1CF 3 of 
the bulge mass, and galaxy M87 is though t to have a supermas- 
sive black hole of mass (3 ± 1) x 10 9 M Q tearing & Rix l2004l) . 
This indicates that by ascertaining a fourth lensed image, a relevant 
upper bound on the mass of the supermassive black hole can be 
obtained. Note also that this method is viable at cosmological dis- 
tances which are inaccessible to precise mass measurements from 
stellar kinematics. 

However, this approach has two main difficulties: first, it re- 
quires knowledge of ko and hence the mass and size of the lens, 
and secondl y, the four th image may be missed or be misinterpreted. 
iBowman et al. I l2004t) have recently considered the lensing proper- 
ties of a cored and a broken power-law isothermal sphere with cen- 
tral point mass. Because of their diverging masses, these models 
are only valid centrally and hence a black hole mass bound from 
lensing, as described here, is not discussed. Based on their numer- 
ical simulation of lensing cross-sections, the authors conclude that 
observations of the fourth image, as characteristic lensing signature 
of a supermassive black hole, should require the next generation of 
radio telescopes but nevertheless recommend searching already. 



3.4 Magnification Invariant 

Using the same conventions as in §2.3, the lens equation I32i can 
be recast thus 



C = z- 



1 + zz 



ao 

z 



such that 



/(*) 



+ 



(46) 



(47) 



Although the central black hole has yielded a further image, it is 
still true that the sum of the signed magnifications is unity. 



4 THE PLUMMER LENS WITH BLACK HOLE AND 
SHEAR 

4.1 Lens Equation 

Circular symmetry is however normally broken in gravitational 
lensing. This may be caused by flattening of the lensing galaxy 
itself or by the presence of an external gravitational field, for in- 
stance, in a galaxy cluster. For simplicity, let us consider a constant 
external mass surface density K ex t and shear 7 ex t- Let the Cartesian 
coordinate axes be in the principal directions so that the perturba- 
tionof the lens equation is in diagonal form ("e.g.. lSchneider et al. I 
Il999t p. 250) 



Ay = 



^ext ~T~ 'T'ext 



Ti o 
o r 2 



o 

^ext '"/ext 



(49) 



where Ti, T2 are the components of the total external shear. The 
lens equation 1321 is now modified thus 



2/i 



1)2 



X 2 



K0X1 



l + ( — ) 

KQX2 
l + ( — ) 



ax 1 



ax 2 
~^2~ 



r 2 x 2 . 



(50) 



(51) 



For Ti 7^ T2, the imaging equation is a fifth order polynomial in 
zz. However, as we will see in the next section, the maximum num- 
ber of images of a Plummer model with black hole and shear is 6. 
Hence in general, there are spurious roots to the imaging polyno- 
mial and a magnification invariant does not exist. 



so again zf — > as \z\ — > 00 and / has no branch cut. The maxi- 
mum number of images is four. Likewise, the imaging equation is 
of fourth degree in z so there are no spurious roots. Therefore, the 
magnification invariant holds for the Plummer model with central 
point mass as well, 



E 



(48) 



4.2 Image and Caustic Properties 

4.2.1 Astroid Caustic 

When Yi = Y2, eqns (50) and (5l) remain circularly symmet- 
ric with three critical circles, two caustic circles and a critical and 
caustic point, as in §3. 

Consider now the case Ti / r 2 , Ti > OAl^ > 0,thatis,the 
external mass surface density /t ex t > 7cxt- Since, by construction, 
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There is now one pair of cusp points on each of the y\ , 1/2 -axes giv- 
ing rise to an astroid caustic. By crossing the astroid inwards, the 
number of images increases from two to four. If the point source 
is at y = 0, an additional fifth image at the centre of L is pro- 
duced, since by l50l and 15T1 x = is a solution. Note that, in this 
case, circular symmetry is broken and an Einstein ring image can- 
not form. The astroid increases in size as Ti , T2 become more dif- 
ferent, eventually intersecting the already existing caustic curves. 
Caustic domains are thereby created in which the image multiplic- 
ity increases from a maximum of four to six. Such a configuration 
of critical and caustic curves is illustrated in the upper and middle 
panel of figure|5] 



4.2.2 Metamorphoses 

As the difference in the shear components increases yet further, 
Ti > A T'2 < and the external mass surface density K cxt < 
7 cx t. Now critical curves and caustics rejoin, as can be seen in the 
lower panel of figure|3| as the inner two critical curves in L merge 
on the 12-axis, the tangent vector vanishes and V det J — 0. This 
corresponds to the rejoining of the inner two caustics on the j/2-axis 
to give two cusps, and is an example of a beak-to-beak metamor- 
phosis. 

Our model can also serve to illustrate the metamorphosis of an urn- 
bilic. Referring to the deflection potential ^ in equation <31> and 
the lens equation with external s hear 1501 . 15 H . the conventional 
form (e.g., Schneider et al. Il999t p. 162) of the Jacobian J of the 
lens mapping is, 



Figure 3. Plummer lens with point mass and external shear, ko = 2, ao = 
0.02. The difference between the components of external shear increases 
from top to bottom. Upper left: critical curves, upper right: caustics, for 
Ti = 0.5, T2 = 0.4. Middle left: critical curves, middle right: caustics, 
for Ti = 0.5, T2 = 0.1. Lower left: critical curves, lower right: caustics, 
for Fi = 0.5, T2 = —0.35. Number labels indicate number of images. 



the coordinate axes point in the principal directions of the exter- 
nal shear, the critical curves in L intersect the si-axis with a tan- 
gent vector parallel to (0, 1) and the X2-axis with (1, 0). The cor- 
responding caustic tangent vanishes if, in addition, dyi /dx2 = 
or dyi/dxi = 0, resp ectively. Then, th ese are cusp points on the 
2/i,2/2-axes in S (e.g. 1 Schneid er et al. Ill999l) . p. 252). They may 
be categorised thus: 



• Cusp points on yi-axis: X2 

K 1 



OA 



= 



= 



2/1 



1 - 

xi(r 2 



r 2 i + (f^) 2 

-Ti) from 



1 - r 2 x\ 



from (15 II 



There is precisely one solution for x\, hence precisely two for xi 
and two for yi since Ti 7^ IV 



• Cusp points on 7/2-axis: x\ 
1 



OA 



9yi 



= 



= 



2/2 



from 1501 1 



x 2 (Ti -T 2 ) from J3TJ 



Again, there is precisely one solution for x%, hence precisely two 
for X2 and two for j/2 since Ti 7^ T2. 

Therefore, Ti 7^ T2 lifts the degeneracy of the caustic point at 
y = which was found the previous circularly symmetric cases. 



J-- 
where 

71 = 

72 = 



1 — K — 7l 
-72 

1 ( dH 

2\dx\' 

d 2 t> 

dxidx 2 
1 ( 9 2 * 
2 



-72 

K + 71 



(52) 



a 2 * 



dx\ ~*~ dx\ 



+ 7cxt 



+ Ke 



The umbilic point occurs where both eigenvalues Ai, A2 of J van- 
ish, 



Ai 
A 2 



1 — K + V71 + 7cxt = 
1 - K - V71 + 7cxt = 0. 



These two conditions can now be applied to derive the location 
of umbilic metamorphoses. Since, by construction, the coordinate 
system in the lens plane is aligned with the principal directions of 
external shear, we can set Ti > T2 and X2 = 0. Then keeping 
Ko, ao, Ti fixed defines a one-parameter family of lens models in 
r 2 . Letting £1 = Xi/xo, we hence obtain the value of T2 for the 
metamorphosis from 



(1+d 2 ) 2 

K 



r 2 



Figure |4| illustrates the caustics for a particular Plummer model 
with central point mass and external shear undergoing metamor- 
phosis at an umbilic point. 
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Figure 4. Caustics for the Plummer lens with central point mass and external shear, kq = 2, ao = 0.06, Ti = 0.6. The difference between the components 
of external shear increases left to right as Y2 decreases. Left: Y2 = —0.1, middle: Y2 = —0.52 where the umbilic metamorphosis occurs, right: Y2 = —0.9. 
Number labels indicate number of images. 



5 CONCLUSION 

This paper provides an analytic study of the lensing properties of 
the Plummer model, together with modifications incorporating the 
effects of a central point mass and external shear. This is a simple 
model for lensing by a flattened galaxy with a central black hole. 
The maximum number of images given by the circular Plummer 
lens is 3. This number rises to 4 if a central point mass is introduced 
and to 6 if external shear is added as well. 

There are three new results of our work. First, a magnification 
invariant is found to exist for both the Plummer lens alone, and for 
the Plummer lens with black hole. Even though the introduction of 
a black hole changes the maximum number of images from 3 to 
4, it is nonetheless true that the sum of the signed magnifications 
remains unity. 

The second result is a treatment of the caustic structure for 
the Plummer lens with black hole and shear. Without shear, there 
are two caustic circles and a caustic point. This gives three disjoint 
regions in the source plane where double, quadruple and double 
imaging respectively occur. The introduction of shear causes the 
caustic point to unfold as an astroid caustic. The presence of four 
cusps on an astroid caustic is established, demonstrating the change 
of image multiplicities as the caustics intersect and metamorphoses 
occur. For modest shear, the astroid is contained within the two 
outer caustics. However, as the components of shear Ti and T2 be- 
come more different, the size of the astroid increases and additional 
caustic domains are created in which the image multiplicity is 6. 

Based on this discussion of caustic curves, we describe a pos- 
sible method to obtain an upper limit on the mass of a black hole 
from gravitational lensing. As this approach relies on counting im- 
ages alone, it could be used at cosmological distance and there- 
fore as an independent check on mass estimates obtained from 
other methods, such as stellar kinematics and reverberation map- 
ping. This emphasises the value of lensing surveys for constraining 
the mass distribution in the inner parts of galaxies. 
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